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ABSTRACT 

Improved photometric sensitivity from space-based telescopes have enabled the detection of phase 
variations for a small sample of hot Jupiters. However, exoplanets in highly eccentric orbits present 
unique opportunities to study the effects of drastically changing incident flux on the upper atmospheres 
of giant planets. Here we expand upon previous studies of phase functions for these planets at 
optical wavelengths by investigating the effects of orbital inclination on the flux ratio as it interacts 
with the other effects induced by orbital eccentricity. We determine optimal orbital inclinations for 
maximum flux ratios and combine these calculations with those of projected separation for application 
to coronagraphic observations. These are applied to several of the known exoplanets which may serve 
as potential targets in current and future coronagraph experiments. 
Subject headings: planetary systems - techniques: photometric 
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1. INTRODUCTION 

The changing phases of an exoplanet as it orbits 
the host star has long been consid ered as a means for 
their detection and chara cterization (jCharbonneau et al.l 
119991 iLeigh et al.1 l2003[ ). The relation between gi- 
ant planet atmospheres and phase curv es have been 
described in d etail by ISudarskv et al.l ((2005) and 
iKane k Gelinol (|2010t ). hereafter KG 10. However, the 
relatively small flux ratio of the planet to the host star 
has presented a major hinderance to the realization of 
such detections. A new era of optical and Infra-Red 
(IR) telescopes are enabling phase detections in a manner 
which was previously inaccessible from the ground. 

Attempts to detect phase signatures have primarily 
been for transiting planets. Examples of observed phase 
variati ons in the IR inclu de HP 189733b (iKnutson et all 
I2009al) and HD 149026b (|Knutson et al.ll2009ri) . Exam- 
ples i n the optical includ e Kepler observations of HAT- 
P-7b C W elsh et al.l l20ldh and phas e variations detecte d 
in the light curve of CoRoT-lb ([Snellen et al.l I2009D . 
Phase variation detection of non-transiting planets have 
been restricted to hot Jupiters, includi ng v And b 
(iHarrington et al.ll2006h and HD 179949b ([Cowan et al.l 
|2007|) . There has also been recent observational evidence 
for the detection of ph ase variations b y the n on-transiting 
planet HD 46375b by iGaulme et all (12010ft. Planets in . 
eccentric orbits, s uch as HD 1715 6b (Barb ieri et al.ll2~007f) 
and HD 80606b (|Laughlin et al.ll2009D . will produce rel- 
atively high phase amplitudes during a brief interval. 

The orbital inclination is usually preferred to be edge- 
on for optimal detection of phase variations since this 
ensures that full phase will occur along the line-of-sight. 
However, this is not necessarily the case depending upon 
the eccentricity and periastron argument of the orbit, 
and indeed there are cases where the maximum flux ra- 
tio occurs when the orbit is face-on. Even though the 
radial velocity technique is biased towards the detection 
of edge-on orbits, since these have a larger radial veloc- 
ity semi- amplitude, certainly this is not alwa ys the case. 
Astrometric studies such as that performed bv lHan et al.l 
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(|2001f) have shown that there are a variety of orbital in- 
clinations as expected, and so it is prudent to consider 
this for attempted phase detections. 

Here we present a thorough exploration of orbital 
parameter-space and the impact upon flux ratios for ec- 
centric planets. The main causes of degeneracy are the 
periastron argument and inclination for a given eccentric- 
ity. We further calculate projected separations at apas- 
tron as a function of inclination and determine their cor- 
respondance with maximum flux ratio locations. With 
all orbital elements considered, one may then use the 
phase locations of maximum flux ratio and maximum 
projected separation to quantify the suitability for follow- 
up observations. 

2. FLUX RATIO COMPONENTS 

In this section, we outline the theoretical frame- 
work which will be used through the remainder of 
th e paper. This formalism has been used extensively 
bv iCollier Cameron et all (j2002[ ) and more recently by 
iRodler et al.l (|2010l) . We refer the reader to KG10 for 
a more detailed description of the particular formalism 
used here. 

The flux ratio of a planet with radius R p to the host 
star is defined as 



e(a, A) 



M^= Ag (\)g(a,\)^- (1) 



where the flux is measured at wavelength A. This flux 
ratio consists of three major components; the geomet- 
ric albedo A g (X), the phase function g(a,X), and the 
inverse-square relation to the star-planet separation r. 
Each of these are briefly described below in the context 
of eccentric orbits. Note that the observed flux ratio from 
an exoplanet is wavelength dependent in that the atmo- 
spheric composition drives the scattering properties and 
thus the shape of the geometric albedo and phase func- 
tion. As was the case with KG10, we confine our study 
to optical wavelengths centered on 550 nm. 

Atmospheric models have shown that there is a de- 
pendence of the geometric albed o of giant planets on 
the semi-major axis of the orbit ( Sudars kTeTaTI [2000, 
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2005; ICahov et al.ll2010D . To account for this, we use the 
analytic function described by KG10 which was i n turn 
derived from the models of ISudarskv et al.l ((2005) . This 
function results in a time-dependence of the albedo as the 
strong irradiation of the atmospheres of giant planets re- 
moves reflective condensates from the upper atmospheres 
during periastron passage. This function can have a dra- 
matic effect in dampening the flux ratio at small star- 
planet separations. 

The phase angle a is defined to be zero when the planet 
is at superior conjunction and is described by 

cos a = sin(a; + /) (2) 

where uj is the argument of periastron and / is the true 
anomaly. Thus, in terms of orbital parameters, minimum 
phase occurs when uj + f — 90° and maximum phase 
occurs when ui + f = 270°. 

The exact nature of a planetary phase function de- 
pends upon the assumptions regarding the scattering 
properties of the atmosphere. Rather than assuming 
isotropic scattering (Lambert spher e) , we adopt t he em- 
pirically derived phase function of iHiltonl ()1992l ) which 
is based upon observations of Jupiter and Venus and 
incorporates substantially more back-scattering due to 
cloud-covering. This approach contains a correction to 
the planetary visual magnitude of the form 

Am(a) = 0.09(a/100°) + 2.39(a/100°) 2 - 0.65(a/100°) 3 

leading to a phase function given by 

g(a) = io-°- 4Am ( Q ) (4) 

which is used th roughout the remainde r of this paper. 

As shown by ISudarskv et~all (|2005[ ) and KG10, the 
maximum flux ratio does not necessarily occur at zero 
phase angle for a non-circular orbit. This is because the 
star-planet separation r is constantly changing and is 
given by 



where a is the semi-major axis and e is the orbital ec- 
centricity. The Rp/r 2 component of Equation [T] becomes 
dominant for highly eccentric orbits which is especially 
important when the orbital inclination is considered. 

3. ORBITAL INCLINATION 

To add the effect of inclination angle to the phase func- 
tion, the phase angle (Equation ^ is modified as follows: 

cos a — sin(w + /) sin i (6) 

At 1st and 3rd quarter {a = 90° and a = 270°), the 
flux ratio is completely independent of inclination angle. 
However, at all other phase angles there arises a complex 
pattern of flux ratios from the inclination and argument 
of periastron. Here we describe these dependencies. 

3.1. Edge- on/ Face- on Orbits 

The two cases worth considering first are the extremes 
of edge-on (i — 90°) and face-on (i = 0°) orbits. As 
noted earlier, the radial velocity method is biased to- 
wards the detection of edge-on orbits due to the increase 
in the semi-amplitude of the signal. In addition, most of 
the planets monitored for phase signatures are known to 



transit. Thus the case of edge-on orbits is currently the 
dominant form of investigated systems. This guarantees 
the observability of both zero phase and full phase, the 
contribution of which to the flux ratio depends upon the 
star-planet separation at a = 0°. 

The case of face-on orbits means that the phase func- 
tion becomes completely flat since only half phase will be 
visible at any one time. Thus the flux ratio is completely 
determined by the eccentricity of the orbit which drives 
both the star-planet separation and the changing albedo 
of the upper atmosphere. 

3.2. Generalized Orientation 

Beyond the cases of edge-on and face-on orbits, the in- 
teraction of inclination and periastron argument becomes 
more complex. Figure [I] demonstrates this for four eccen- 
tric orbits, each with fixed periastron arguments, which 
are inclined from edge-on to face-on. Within the range of 
45° < uj < 135° the flux ratio actually increases with de- 
creasing orbital inclination as the dayside of the planet 
becomes more visible where the planet is the hottest, 
most noticable for the case of uj — 90° . Outside of this 
regime, the peak flux from the planet generally declines 
as the inclination increases and thus the access to the full 
phase diminishes. 

3.3. Peak Flux Ratio Maps 

A further level of detail to the description of inclination 
dependence may be added by calculating the peak flux 
ratio for the full range of inclinations (0° < i < 90°) and 
periastron arguments (0° < uj < 360°) for a given period 
and eccentricity. The resulting intensity maps then show 
the optimal orbital configuration for detection and how 
the distribution of peak flux ratios smoothly varies with 
these configurations. 

Shown in Figure [2] are two such examples of these in- 
tensity maps for eccentricities of 0.2 and 0.8. In each 
case, the strongest flux ratio occurs where the full phase 
of the planet coincides with the smallest star-planet sep- 
aration (i = 90°, uj = 270°). However, notice the in- 
teraction which occurs where uj — 90°. The changing 
peak planetary flux along this line of the intensity maps 
is due to the competing components of the phase func- 
tion and star-planet separation as the planet moves from 
a crescent phase near periastron to the full phase of the 
planet at apastron. For highly eccentric orbits, the planet 
flux at uj = 90° becomes the brightest when the orbit is 
viewed face-on since the flux ratio is dominated by the 
star-planet separation. 

4. PROJECTED SEPARATION 

The star-planet projected separation is a component 
which will influence the target selection for coronagraph 
experiments for direct detec t ion o f the reflected plan- 
etary flux. iBeichman et al.l (|2010t ) tabulate the inner 
working angle of selected future ground-based imaging 
instruments which range from 0.03" to 0.17", compared 
to 0.035" to 0.850" for James Webb Space Telescope 
(JWST) instruments. The angular projected separation 
of the planet from the star is given by 

A6 = ~ (cos 2 (a; + /) + sin 2 (w + /) cos 2 i) * (7) 

where d is the star-observer distance. Here we utilize the 
analysis of the previous section to determine the orbital 
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Fig. 2. — Intensity maps for a P = 100 days planet with e = 0.2 (left) and e = 0.8 (right). These maps show the peak flux ratio values 
for the full range of inclinations and periastron arguments. For e = 0.2, the minimum peak flux ratio occurs where u) = 90° and i = 13°. 
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location of maximum flux ratio for a subset of the known 
exoplanets. In Table [TJ we report the true anomaly / 
(angle between the direction of periapsis and the current 
position of the planet in the orbit), flux ratio e, and pro- 
jected separation A9 at this location for the most eccen- 
tric exoplanets. We additionally report if the maximum 
flux ratio occurs for an edge-on or face-on orientation. 
Notice that for face-on orbits the true anomaly is zero at 
this location since the flux is completely driven by the 
event of periastron passage. 

There is a bias against optimal separation since maxi- 
mum flux naturally occurs when the planet is closest to 
the star. This is most severe for edge-on orbits, such 
as HD 80606b where the angular separation is ~ 0" 
at this location. For this reason, the location of max- 
imum flux ratio preferentially corresponds to the loca- 
tion of minimum angular separation. Consider the case 
of HD 39091b (Figured]) where the peak flux arises from 
an edge-on orientation. The peak flux ratio corresponds 
to a 0.02" separation, whereas the planet reaches a sep- 
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HD 39091 b 
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Orbital Phase 

Fig. 3. — The flux-ratio (solid line) and star-planet angular 
separation (dashed line) for HD 39091b, assuming an edge-on orbit. 

aration of ~ 0.26" during the entire phase. By matching 
the angular separations above the resolution criteria of 
the coronagraph, one can choose the optimal targets for 
monitoring during predicted the orbital phase of peak 
flux ratio. 

5. EXOPLANET CHARACTERIZATION 

The main targets for which phase signature detection 
will be attempted will likely be those for which a planet 
is already known to be present. In this case, the pur- 
pose of the observations are chatacterization rather than 
discovery. As shown by Equation [TJ the flux ratio is a 
function of the geometric albedo A gi the periastron ar- 
gument a/, the orbital inclination i, and the planetary 
radius R p . The challenge of disentangling those com- 
ponents from the phase shape and amplitude will vary 
depending upon what is already known for that system. 

If the planet is known to transit then one can measure 
i, uj, and R p . One can thus determine A g which, along 
with models of the planetary structure based upon the 
size and mass, can constrain the properties of the at- 
mosphere. If however the planet does not transit, then 
we can estimate the radius of the planet from the mea- 
sured mass using theoret ical models, such as tho se of 
iBodenheimer et all (|2003t) and lFortnev et al.l (|2007f ). for 
which there will be associated uncertainties depending 
upon factors such as the age of the planet and the as- 
sumed core model. In some cases, astrometric measure- 
ments of t he host star can p erformed, such as those car- 
ried out bv lHan et al.1 (|200 lj). This will have the simulta- 
neous result of resolving the inclination of the orbit, thus 
breaking the degeneracy of the flux ratio with the peri- 
astron argument, and determining the true mass of the 
secondary companion. For non-transiting planets discov- 
ered using the radial velocity method, one has the advan- 
tage of a strong bias towards bright host stars relative to 
those discovered using the transit method. The increase 
in signal-to-noise will be a great asset for an appropriate 
instrument to exploit, as described in Section [6] 

Possible contamination of the phase signature at pho- 
tometric passbands by internal heating of the planet is 
worth considering. If indeed the thermal signature of 
the planet is dominated by internal heating rather than 
incident flux from the host star, the signature will main- 
tain a constant offset of the predicted photometric sig- 
nature. The self-luminous properties of the planet will 
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be a strong function of the age of the system, and will 
be restricted to planets whose age is (< 1 Gyr). Cur- 
rent radial velocity surveys choose targets based upon 
stability as well as spectral line features and thus pre- 
fer F-G-K ma in sequence st ars where activity is known 
to be reduced (|Wrightl 120051 . Thus, targets drawn from 
the radial velocity discoveries will preferentially be of an 
age whereby thermal contamination of the optical phase 
signature will be minimal. 

6. FEASIBILITY REQUIREMENTS 

In this section we briefly outline requirements which 
will influence the detection of the previously described 
phase signatures. 

6.1. Instrumentation 

The detection of the phase amplitudes described here 
present a significant challenge to instrumentation re- 
quirements. As mentioned in Section [1] optical phase 
variations have been detected in Kepler photometry. 
However, the vast majority of Kepler targets are rela- 
tively faint and so not ideally suited towards character- 
ization with current follow-up capabilities. Conversely, 
planets discovered through radial velocities have bright 
host stars by comparison. 

The instrumentation requirement for successful detec- 
tion of the signatures shown in Table [1] is photometry 
with an accuracy of ~ 10 -6 . The necessary stability of 
this precision over long timescales depends upon the frac- 
tion of the orbital phase over which the largest change in 
phase amplitude occurs, described in detail by KG10. 
The photometer for the Kepler mission, for example, 
is designed to achieve high-precision photometry over 
the 6.5 hour window of a transit, but is not designed 
for long-term stabili ty over the lifetime of the mission 
(|Borucki et al.ll2~010l ). Fortunately the orbits of the ra- 
dial velocity planets are well understood in most cases 
and so we can accurately predict not just the amplitude 
of the predicted phase signature but also the phase and 
times of maximum and minimum flux ratios. This knowl- 
edge will help to distinguish the phase signatures from 
instrumental drift effects. 

From the ground, the challanges are more substantial 
since one needs to also contend with the offsets from 
night-to-night variations. Future generation telescopes 
will provide opportunities to achieve very high preci- 
sion, such as the European Extremely Large Telescope 
(E-ELT), the Thirty Meter Telescope (TMT), and the 
Giant Ma gellan Telescope (G MT). It has been demon- 
strated bv lColon etail ()2010D that precision photometry 
of < 0.05% can be achieved with large telescopes through 
the use of narrow-band filters. 

The phase variation due to the planetary orbit is the 
main observable we have discussed here but these cal- 
culations can be combined with projected separations 
for use with coronagraphs, as described in Section 01 
Ground-based coronagr aph experiments, suc h as that be- 
ing pursued at Palomar (Hinkley ct al. 20ll|), are rapidly 
improving in instrument efficiency and the suppression 
of speckle noise. Simulations of these coronagraphs in- 
dicate that long-term stability may indeed be possible 
(|Beichman et al.ll2010f) . though this may only be feasi- 
ble in the short-term at the needed precision for young 
self-luminous planets. Space-based instruments, such as 



the planned Lyot coronagraph on NIRCam for the James 
Webb Space Telescope (JWST), may be able to achieve 
phase detections for a sample of the most favorable tar- 
gets, though in this case the instrument is optimized to- 
wards young planets around late-type stars. 

6.2. Stellar Variability 

At the level of photometric precision required here, 
it is important to consider the level of intri nsic stellar 
variab ility. An analysis of Kepler data by ICiardi et al.l 
(2010) found that most dwarf stars are stable down to 
the the precision of the Kepler spacecraft, with G-dwarfs 
being the most stable of the studied spectral types. The 
well-known orbital parameters of the radial velocity tar- 
get stars will aid in separating the signals of planetary 
phase from that of the host star variability. Addition- 
ally, most of the known exoplanet host stars have been 
well-characterized through extensive ground-based pho- 
tometry and spectroscopy. 

The main cause of photometric variability in F-G-K 
stars is starspots and rotation, as verified by t he Ke - 
pler variability study performed by Bas ri et all (|2011h . 
The effects of starspots on exoplanet dete ction have been 
previo usly studied, such as the work of iMakarov et all 
(2009). Since the distribution of the starspot periodicity 
is related to the rotation rates of these stars, the detec- 
tion of short-period planets (such as most of the known 
transiting planets) is relatively unaffected by these vari- 
ations. However, there will inevitably be cases where the 
orbital period of the planet is close to the rotation period 
of the star, which is gener ally in the range 10-4 days for 
radial velocity host stars (jSimpson et al.|[2010l ). In such 
cases, the peaks in the power spectrum from a fourier 
analysis of the photometry may separate to a degree 
where the starspot variability can be isolated from the 
phase signature. It should be noted that disentangling 
these signals may substantially increase the required ob- 
serving time. 

7. CONCLUSIONS 

The push towards characterizing the atmospheres of 
exoplanets will rapidly expand once further direct mea- 
surements of atmospheric albedos and thermal properties 
become possible. Current space telescopes are already 
detecting exoplanet phase variations in the optical (eg., 
Kepler) and the IR (eg., Spitzer). One approach to opti- 
mizing searches for phase variations in known eccentric 
systems is the refinement of or bital parameters t hrough 
radial velocity measurements (|Kane et~all[2009l) . The 
peak flux ratio maps may then be used to constrain the 
detectability of each system and the location of the peak 
intensity with respect to the projected separation, par- 
ticularly important for proposed coronagraph missions 
such as the Spectro-Polarimetric Imaging and Charac- 
terization of Exo-planetary Systems (SPICES) mission. 
More imminently, Kepler will soon detect transiting long- 
period (P > 100 days) planets where the bias will be 
towards eccentric or bits since those have a hig her prob- 
ability of transiting (jKane fc von Braunll2008f ). In this 
case the results of the presented work may be utilized to 
predict flux ratios for edge-on orbits as an independent 
measurement of the periastron argument. The faintness 
of these host stars will provide many instrumentation 
and observing challenges which make known radial veloc- 
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ity planets more attractive targets despite their unknown 
orbital inclination. Even so, as more science results are 
released by the Kepler mission, the study of photometric 
phase variations of long-period planets will become an 
increasingly relevant endeavour in the characterization 
of exoplanets. 
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